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Abstract. Hambly, Keevash, O'Connell and Stark have proven a central limit 
theorem for the characteristic polynomial of a permutation matrix with respect 
to the uniform measure on the symmetric group. We generalize this result in 
several ways. Wc prove here a central limit theorem for multiplicative class 
functions on symmetric group with respect to the Ewens measure and compute 
the covariance of the real and the imaginary part in the limit. We also estimate 
the rate of convergence with the Wasserstein distance. 
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1. Introduction 

The study of random matrices has gained importance in many areas of mathematics 
and physics, for example in nuclear physics, infinite dimensional integrable systems 
and large— n representation theory. Random matrix theory (RMT) was in the recent 
years also of big interest in number theory since the study of the spectrum of the 
characteristic polynomial of a random matrix in a compact Lie group was central 
in obtaining conjectures. A good example to illustrate this is the paper of Keating 



and Snaith 1 10 . They conjectured that the Riemann zeta function on the critical 
line could be modeled by the characteristic polynomial of a random unitary matrix 
considered on the unit circle. One of the results in ||l^ is 

Theorem 1.1. Let x be a fixed complex number with |a;| = 1 and g„ be a unitary 
matrix chosen at random with respect to the Haar measure. Then 



Log dct(/„ - X5„) 

^ ■ Afi + iAf2 for n ^ oo (1.1) 



i log(n) 

and Ml , M2 independent, normal distributed random variables. 

Constin and Lebowitz have proven 10 years earlier in |6| a weaker version of this 
theorem. They showed that 



'Log(^det(/„ - xgn)j 



Im( \ ^ I 4M. (1.2) 

i log(n) 



They conjectured that the same is true for the real part and that the imaginary 
part and the real part are independent in the limit, but haven't been able to prove 
this. 

The situation for the characteristic polynomial of permutation matrices is similar. 
A permutation matrix is a unitary matrix of the form ((5i.o-(j))i<i,j<n with a € Sn 
and Sn the symmetric group. It is easy to see that the permutation matrices form 
a group isomorphic to Sn- We call for simplicity both groups Sn and use this 
identification without mentioning it explicitly. The characteristic polynomial Zn {x) 
of a permutation matrix is defined as 

Znix) = Zn{x){a) := dct(/ — xa) with a; G C, cr £ Sn- (1-3) 

Hambly, Keevash, O'Connell and Stark have proven in Q 

Theorem 1.2 (B.M. Hambly, P.Keevash, N. O'Connell and D.Stark). Let g G Sn be 

chosen uniformly at random and x be a fixed complex number with \x\ — I, not a 
root of unity and of finite type (see Definition 3.1^ . Then 



Log(Z„(a;)) \ d , / Log(Z„(a;)) 



Re I "7^"'"^^^ I AAA,, Im I "7^"""-^^ I 4 A4 (1.4) 
/f|log(n)/ \^^\og{n) 

with Ma , Ab standard normal distributed random variables. 

As for unitary matrices, it is natural to ask if Ma and At are independent. This 
question is not considered in and was the main motivation for this paper. We 



will see in Corollary 2.7.1 that Ma and Mb are indeed independent. But we can 



show here much more. The main result of this paper is Theorem 2.7, which is an 



extension of Theorem 1.2 with three important differences. These differences are 
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We compute the covariance of the real and the imaginary part in the hmit. 



We endow S'„ with the Ewens measure (see Definition 2.4), which is a gen- 
erahzation of the uniform measure. 

We consider more general class functions on Sn, the so called multiplicative 



class functions W"{f) (see Definition 2.5) 



We have introduced multiplicative class functions M^"(/) in as generalization of 
Zn{x) and studied there the asymptotic behavior of their moments with respect to 
the uniform measure on Sn- 

The main idea there is to write down the generating function with a combinatorial 
argument and to use function theory to extract the asymptotic behavior. 

The structure of this paper is as follows: we introduce in Section |^ multiplicative 
class functions and state the main theorem of this paper. We do in Section ^ 

some preparations and state in Section ^ an auxiliary central limit theorem. We 
then prove in Section ^ the main Theorem 2.7. In Section |^ we then estimate 



the convergence rate of the probability measures in the main theorem with the 
Wasserstein distance. 



2. Definition and main theorem 



We introduce in Section 2.1 the Ewens measure and some well know functions on Sn- 
In Section 2.2, we give an alternative expression for Zn{x) and use this expression to 



introduce the multiplicative class functions W"{f){x). We then state in Section 2.3 
the main theorem of this paper. 



2.1. The symmetric group Sn- All functions on Sn in this paper are invariant 
under conjugation {u{hgh~^) = u{g)), i.e. they are class functions. It is therefore 
natural to take a look at the conjugation classes of Sn- These can be parameterize 
with partitions. 

Definition 2.1. A partition X is a sequence of nonnegative integers Ai > A2 > ■ ■ ■ . 
The size of the partition is |A| := Am dnd the length of £{X) is the largest £ such 
that Xi 7^ 0. We call X a partition of n if \X\ = n, and denote this by X\- n. 

Let a (z Sn he arbitrary and write a = aia2 ■ • ■ cfi € Sn with ai disjoint cycles of 
length Ai. Since disjoint cycles commute, we can assume that Ai > A2 > • • • > A^ . 
We call the partition A ~ (Ai, A2, . . . , A^) the cycle type of a. It remains to show 
that two elements of 5„ are conjugated if and only if they have the same cycle type. 
Since this is well known, we omit the proof here and refer to chapter 39] . 
We introduce as next the cycles counts of a given length. 

Definition 2.2. Let a ^ Sn be given with cycle-type A ~ (Ai, • • • , A^). We define 
Cm C'^' (a) := # {^ : 1 < z < / and A, = m} . (2.1) 

The functions Cm ^ [a) depends only on the cycle-type of a and are thus class func- 
tions on Sn- It is clear that the cycle type A of tr S S'„ is uniquely determined by 

the values of C^"-* , • • • , Ci""* . We therefore can work with the functions Cm ' or with 

ill) 

the cycle type A. We prefer here to use Cm - 

The most natural measure on Sn is the uniform measure (i.e P [A] = -1^), but 
there exists of course other important measures on Sn than the uniform measure. 
One of these measures is the Ewens measure, appearing in population genetics (see 
1^). The Ewens measure is a generalization of the uniform measure and has an 
additional weight depending on the total number of cycles. 
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Definition 2.3. Let 9 > 0. We then set for a e Sn with cycle type A 

' '-^ e{e + i)---{n + e-i)- 

The measure ¥g [.] is called the Ewens measure with parameter 6. 



(2.2) 



The uniform measure is the special case 6 = 1. 

Many resuhs about the Ewens measure can be found in the book chapter 4] , for 
instance 

Lemma 2.4. Let 6 > be given. The distribution of (Ci, C2, ■ ■ ■ , C„) with respect 
to the Ewens measure on Sn is given by 

P.[(Ci = ci,...,C„ = c„)] = -^^n— (-) lU=E'^^™ (2-3) 
and the expectation of cln' is 

„ (9+n-m-l 



E 



^ \{m<n). (2.4) 



n—m 



2.2. Definition of VK"(/). We now come multiplicative class fonctions. We first 
give a explicit expression for Zn{x) using the cycle counts Cm \ and then use this 
expression to introduce multiplicative class functions. 

Let a ^ Sn he given, then 

n 

Zn{x)=l[{l~x^f'r:^\ (2.5) 

m— 1 

The proof of this equation is straightforward. One first has to take a look a the 
case when cr is a cycle and then take a look at the general case. More details can 
be found in p^ . 

Definition 2.5. Let x E C be complex number and f : S^ ^ C be a real ana- 
lytic function with S^ := {z G C : \z\ = 1}. We then define the multiplicative class 
function associated to the function f as 

n 

W"{f) = W"{f){x){<7) := H /(a;™)^^"' . (2.6) 

7n—l 

For brevity, we simply call this a multiplicative class function. We also set 

n 

w^{f){x) = \og{W^\f){x)) := J2 Ci"^ Mfi^n) (2.7) 

m— 1 

with log the principal branch of logarithm and log(— ?/) := ^og{y) + i-K for y G R>o 
and log(O) = oo. 

It is clear that VF"(/) and Wn{f) are class functions on Sn and that the characteristic 
polynomial is the special case f{x) = 1 — a;. 

Multiplicative class functions W"{f) have been introduced in [Q as generalization 
of the characteristic polynomial Zn{x). The motivation was that one can compute 
the asymptotic behaviour of the moments of Zn{x) and of VF"(/) using the same 
method with only minor changes. We will see in Section ^ that this is here also the 
case. 

In contrast, the extension to the Ewens measure is much more laborious. We will 
see in the proof of Lemma 4.2 that the case < < 1 causes much more work than 
the uniform measure. 
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2.3. The main theorem. We state in this section the main resuh of this paper. 
For this we need some smaU preparations. 

Definition 2.6. Let / : 5^ — > C &e a real analytic function, x £ be arbitrary but 
fixed. If X is a root of unity of order p, i.e. x^ = 1 and p minimal, we set 

"^(/)(^)-^El°g(/(^"))- (2-8) 

P m=l 

// X is not a root of unity, we set 

m{f){x)= f\og{f{e'^n) ds. (2.9) 



// the sum in (2.8) and the integral in ( ^.9| ) respectively does not exists, we set 
™(/)(a;) := oo. 

The integral in (^^) exists for each real analytic / 7^ since the zeros of / are 
isolated and 

log(/(e2--)) ~/vlog(s-so) for s ^ (2.10) 

for So a- zero of /. Thus m{f){x) 00 if and only if / = 0, or a; is a root of unity 
of order p and /(x™) = for some 1 < m < p. 

One could rewrite the integral in (2.9) as J^i log(/((p)) d(p with d(p the uniform 
measure on . We have not used this because this can be easily confused with the 
complex integral J^, log(/(z)) dz = log(/(e2"«))e2"" ds. 

We now come to the main theorem. We distinguish the cases x a root of unit and 
X not a root of unity. For x not a root of unity we assume, as in Theorem |l.2| , that 
X of finite type (see Definition ^3.15 ). This condition is essential in our proof. We 
postpone the definition of finite type to the end of Section 3.3 since we can illustrate 
there why we need this assumption. 

Theorem 2.7. Let f : C be real analytic, x G and Sn be endowed with 

Ewens measure. 

• If X is not a root of unity and of finite type, and all zeros of f are roots of 
the unity, then 

"^^D^ 0VbiRm(/)(a;) 4 N (2.11) 

with M a complex normal distributed random variable. The covariance of 
the real and the imaginary part of M is given by 



— Im 
2 



l^j\og\f{e^-n) ds^ (2.12) 

The real and the imaginary part of M are independent if and only if the 
covariance is equal to 0. 

If X is a root of unity of order p and /(.t™) 7^ for all 1 < m < p then 

^^^5= - 0V^)rnif ) 4 Af (2.13) 
^log(n) 

with Af a complex normal distributed random variable. The covariance of 
the real and the imaginary part of M is given by 

^I-(^El°g'(/(-'")))- (2-14) 
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The real and the imaginary part of M are independent if and only if the 
covariance is equal to 0. 

As promised in the beginning, we now can prove 



Corollary 2.7.1. The random variables Afa and Mb in Theorem l.i are indepen- 
dent. 



Proof. We know that Zn{x) = W^{1 — x) and can therefore can apply Theorem 2.7 



The independence thus follows if we can show that the expression in (2.12) is for 



f{x) ~ 1~ X. A simple computation shows log(l — e = log |1 — e | — ins for 
sG We get 

ilm ^ log'(l - e'""^) d^= j ' Re(log(l - e2"^))lm(log(l - e^^^*")) ds 

.1/2 

= / -sTTlogll-e^'^'^l ds = 0. 
The last integral is since the integrand is odd. □ 

3. Preliminaries 

We present in this section a collection of well known definitions and results, which 



we need for the proof of Theorem 2.7 



in) 

3.1. Asymptotic behavior of C,„ and the Feller-coupHng. We follow in this 
section the book 

Lemma 3.1. Let Sn be endowed with Ewens measure with parameter 9 > 0. For 

(n) 

each m G N, the random variables Cm converge as n —> oo in distribution to a 
Poisson distributed random variable Ym = Ym.e with E \Ym] ~ —. In fact, we have 
for allbGN 

A(yi,r2,--- ,n) (r^^oo), (3.1) 

with all Ym independent. 

Proof See H Section 4]. □ 

One of the problems of convergence in distribution of a sequence (X„)„gN is that 
usually all X„ are defined on different probability spaces and it is thus very difficult 
to compare them directly. This is the case for Cm-* and Cm'^^\ Fortunately, the 
Feller coupling constructs for each 9 > a probability space and new random 

(n) (n) 

variables Cm and Ym on this space, which have the same distributions as the Cm 
and Ym above and can easily be compared. 

The construction works as follows: Let ^ := (1^2C3^4C5 ' ' ' ) be a sequence of inde- 
pendent Bernoulli random variables with E [^m] = g-^^m-i) • m— spacing is a 
sequence of m — 1 consecutive zeroes in ^ or its truncations: 

1 1. 

m— 1 times 

Definition 3.2. Let cln' (^) be the number of m-spacings in 1^2 ' ' ' ^nl- We define 
Ym{C) I'O be the number of m-spacings in the whole sequence ^. 

Theorem 3.3. We have 

• The above- constructed Cm\£,) have the same distribution as the Cm '(A) in 



Definition 



9. 
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• ^Tn(0 finite and Poisson distributed with IK [Y„i{£^)] 

• AllYm{£,) are independent. 



We have 



8(8+1) 



*/^> 1, 

«x i/0<6i<l. 



< J 8+n 
- \ 8(8 + 1 



(3.2) 



• For any fixed 6 € N, 

P[(C(")(0,--- ,C^"HC))^(n(Or-- ,n(0)] ^0(n^cx3). 

Proo/. See |, Chapter 4] and |, Theorem 2]. □ 

We write Cm ^ and for both sets of random variables and do not distinguish 
them anymore. 



3.2. Elementary analysis. We give here some simple and well known results from 
analysis. We state them without further comments. 

Lemma 3.4 (Abel's partial summation) . Let ai , • • • , a„ , 6i , • • • , be complex num- 
bers. Then 

n n — 1 

arnbm ^ri&ii + ^ ^,„(5,„+l - 5„) (3.3) 
m—1 m—1 

wii/i Am := X^fcLi ^/i: 'Ti > 1 and Aq :— 0. 
We can use this to prove 

Lemma 3.5. Let (am)m=i a complex sequence. If 

1 " 

-y am^E + Oin-^) (3.4) 

m—1 

for some S > 0, then there exists a X G R such that 

n 

J2^^E\og{n) + K + 0{n-^), (3.5) 



m 

m—1 



V ^=i?log(2) + 0(n-^). (3.6) 
^ — ' m 

n/2<m<n 



Proof. The lemma follows from Lemma 3.3 and straightforward verification. □ 



Finally we need 

Lemma 3.6 (Holder inequality). Let {am)m=i '^^'^ (^m)m=i be finite sequences and 
Pi 9 ^ 1 such that ^ + i = 1. Then 

n I n \ 1/P / „ \ 1/9 
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3.3. Uniformly distributed sequences. Wc introduce in this section uniformly 
distributed sequences and state some interesting properties. We follow the book 
and omit most of the proofs since they are not difficult. 

We begin with the definition of uniformly distributed sequences. 
Definition 3.7. Let t = (im)m^i be a sequence in the interval [0, 1]. We set 

A„i[a,(3])^An{[a,l3],t) := # {1 < m < n : t,„ e [a,/3]}. (3.8) 
The sequence t — (im)„gjs[ is called uniformly distributed in [0, 1] if 

lim 

n— ^oc 

for each a, (3 with < a < /3 < 1. 

The following theorem shows that the name uniformly distributed is well chosen. 

Theorem 3.8. Let h : [0, 1] — > C 6e a proper Riemann-integrable function and 
t = (i)ri)m=i uniformly distributed sequence, then 



An{[a,P]) 



(3.9) 



- V h{tm) / h{s) ds. 



(3.10) 



Proof. If h{t) = lL[Q,_^](t), then ( ^.10 ) follows immediately from (3.9). For h an 
arbitrary, proper Riemann-integrable function, one can uses an approximation ar- 
gument. □ 



We would like to emphasize at that this point, that ( 3.10 ) does not have to be true 
for improper Riemann integrable functions. One thus needs further assumptions 
to handle functions like log(t) or This is indeed the cause why we need in 

Theorem 2.7 the finite type condition. 

We also need the discrepancy of a sequence. 



Definition 3.9. Let itm)m=i be a sequence in [0, 1]. We then define 



D„ = £)„(t) := sup 

0<a<^J<l 



K = Kit) 



sup 

o<;g<i 



n 

A.([o,/3]; 



(13 -a) 



(3.11) 
(3.12) 



We call Dji the discrepancy and D* the ^-discrepancy of the sequence t. 



It is easy to see that _D* < D„ < 2Z?* and therefore Z?„ and 13* are equivalent. We 
prefer to work with Z3* since we have a more explicit expression for it. 

Lemma 3.10. Let n be fixed and t = (tm),^^^ be a sequence in [0, 1] with ti < 
■ ■ ■ < tr,. Then 



Kit) 



max max 

l<?7i<n 



m 
n 



tri 



1 



(3.13) 



An important fact is that Theorem 3.8, the discrepancy and uniformly distributed 
sequences are closely related. We have 

Lemma 3.11. Let t = (im)m=i be a sequence in [0,1]. The following conditions 
are equivalent 

(1) t is uniformly distributed, 

(2) lim„^ooi?„(t) = 0, 
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(3) Let /i : [0, 1] — > C be a proper Riemann-integrable function. Then 

for n — > oo. 



1 " 

- V h(tr,i) / h{s) ds 



We have introduced t he discrepancy since it aUows us to estimate the rate of con- 
vergence in Theorem |3.8[ We have 

Theorem 3.12 (Koksma's inequahty). Let h : [0,1] ^ C be a function of bounded 
variation V{h), and t — (im),„gN sequence in [0, 1]. Then 



1 " 

- h(tr,i) - / h{s) ds 

ra=\ 



(3.14) 



We will work in Section || with h{s) = Iog(/(e^'^*'') for / real analytic. If sq is a 
zero of /, then 

I I d K 

h(s) ^ ii'log s — So a-nd —h(s) ^ r 

' ' ds \s ~ Sq\ 

for s — >■ So and K a constant. We thus cannot apply Koksma's inequality in this 
situation. We instead use 

Theorem 3.13. Let 6 > and = sq < si < • • • < s^+i = 1 with Sk+i — Sk > 26 
be given. We set 

d 

I -=\Jlk with h ■■= [sk + S, sk+i - 5]. (3.15) 

Further let t = (tm)m=i be a sequence in I and /i : / — > C 6e proper Riemann 
integrable function of bounded variation V{h). Then 



1 " f 
— / h(t,n) — / h{s) ds 
n ^-^ I J 

m=l •' ^ 



< Dl{t)V{h) + 6 \hisk + 6)\ + \hisk+i - 5)1 j 

(3.16) 

with V{h) calculated in I and D* calculated in [0, 1]. 



Proof. We consider here only / = [(5, 1 — 5]. The general case is complete similar. 
W.l.o.g we can assume that ti < t2 < ■ ■ ■ < tn- We put ^o := S,tn+i := 1 — 6 and 
look at 



m=0 * 



s ) dh{s) 



(3.17) 



It follows immediately from Lemma 3.10 that |J| < D^V{h). This observation 
together with a partial integration proves the theorem. □ 



We consider in Section ^ 

1 " 

En := - ^ log(/(e2"™*) as n ^ oo 



(3.18) 



m— 1 



for a fixed t and / real analytic. We are thus interested in the sequence ({n^i})m=i 
with 

{t} t - [t] and [t] := max {n eZ,n<t}. (3.19) 
It is clear that we have to distinguish between t rational and t irrational. The case t 



rational is easy to handle. The next lemma shows that we can apply Theorem 3.12 
and Theorem 3.13 respectively for t irrational. 



10 



DIRK ZEINDLERi 



Lemma 3.14. Lett e K &e given. The sequence ({"^0)m=i uniformly distributed 
in [0, 1] if and only if t is irrational. 

We can use Theorem 3.13| only if we can estimate of the discrepancy of 
and find a S = S{n) with | {mt} — Sk\ > S for 1 < m < n,l < k < d and si, . . . Sd 
the zeros of /(e^'^'*) such that error terms in ( 3.16| ) vanishes for n -> oo. This is 
not so easy to do for arbitrary t and /, but it can be done if t is of finite type and 
ah zeros of / arc roots of unity. 

Definition 3.15. Let x = e^'^'* be given. We call x and t respectively of finite type 
if there exist constants K,^ > Q such that 

\\{mt)\\> ^ for allm(^-L\{0} . (3.20) 

w)«i/i j|sj| := min {s, 1 — s}. We set 

7] = inf {7 G IR+ : 7 fulfills (^^ | . (3.21) 
The constant rj is called the type of x and t respectively. 

A simple computation now shows 

Lemma 3.16. Let t be of type rj and q G N &e fixed. For j > there exist a 
constant K such that 



max 

0<P<<3 



1 n P 
{mt\ 



> — for all m e Z \ {0} . (3.22) 



If t is of finite type, then the discrepancy of the sequence ({™0)m6N '^^^ esti- 
mated with Erdos-Turan-Koksma inequality. We get 

Theorem 3.17. Let t be of finite type rj and t ~ ({'Tit})^^]^ . We then have for 
each e > 



Dr,{t)=0[n~-+^j (3.23) 
Proof. See [0- □ 
4. Auxiliary central limit theorem 

We prove in this section the following auxiliary central limit theorem 

Theorem 4.1. Let 6 > be fixed. Let {cm)m=i ^ sequence of complex numbers 
with a,„ = Re(cm), b,n = Ini(cm) and 

(1) \b,n\ < 27r and |a„J = 0(log(m)), 

(2) ^ Sm=i = Ea + O^n-^") for CO and some 5a > 0, 

(3) ^ Sm=i I''™! = Eb + 0{n~^>') for CO and some Sb > 0, 

(4) ^Em^lflm^K, 7^T.l=lbl^^Vb, ^ E™=1 Om^m £^a6 /or n ^ OO, 

(5) ^j:i=M'=o{loe/'in)), 

(6) It exists ap>\ with i E„/2<m<« 1^™!^ = ^(l)- 

We define A„ = MO) E™=i amC*™^ and B„ = S„(0) := 2™=! b^da\ 
We then have 

An + iBn - Eg [An + iB^] (41) 
Vlog(n) 

with Af a complex normal distributed random variable with covariance matrix 
The real and the imaginary part of N are independent if and only if Eab = 0. 
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Theorem 4.1 is similar to Lemma 3.1 in but there are two important differences: 

• We can calculate the covariance between the real and imaginary part and 
show when they are independent in the limit. 

• We consider a more general measure on Sn- 



Proof. We use in this proof the Feller-coupling (see Section p.l[ ). The random 
variables Cm , Cm and arc therefore defined on the same space. The idea of 

(n) 

the proof is to replace Cm by and to do the computations with Ym- For this 
reason we set 



a,nYm and S„ := ^ b,nY,r, 



(4.3) 



m=l 



We then can show 

Lemma 4.2. Let am,bjn, An, Bn, An and _B„ be as above. We then have 



E 



An + "iBn An iB„ 



0(1) as n -> oo. 



In particular, we see that 
1 



■.{{An+tBn)~{An+tBn)) 



(4.4) 



(4.5) 



and thus . ^ (An + iBn) and , ^ [An + iBn) have the same asymptotic be- 

■y/log(n) V'°S(") 

havior as n oo. 

We first finish the proof of Theorem and then prove Lemma [4.2| . 
We do this by computing the characteristic function of 



.\ (An + iBn). We set 



CXp ita- 



An - E 



An 



V^og{n) 



itb- 



Bn-E 



Bn 



\/log(n) 



(4.6) 



We can compute x^^H'^-a^'^b) explicitly since E [e**^] = exp(A(e'* — 1)) for Y a 
Poisson distributed random variable with E [y] = A. We get 



X^"^ = E 



exp 



itbBn 



= E 




\/^og{n) Vlog(n) 



TO 




An 



itbE 



Bn 



itaE 



An 



itbE 



Bn 



V^og{n) y/log{n) 



\/^og{n) 



E 



(4.7) 



We now use that 



{e^'-l)-{^t-'-J 
1 



1 



1 



< \t^\ and get 



■ {itaa,n + itbbn 



log(n) \ 2'*'"^" ~^ 2^™^'' ~^ ^mbmtatb 



log'^'(n) 



0{\c 



(4.8) 
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for \ta\, \th\ < K with K an arbitrary fixed number. We get 



X^"^(ia,4) =exp 



21og(n) 



E 



log^/2(n) 



-O 



( ^ 1^ 

I log(n) ^ 



m=l 



TO 



(4.9) 



We apply Lemma 3.5 to each summand. The first sunimand converge by condi- 
tion (Q) to —di^t'i + ^tl). The second summand converge by condition (Q) to 
OtatbEab- The third summand converge by the conditions (|), (0) and (||) to 0. 
Therefore 



xn{ta,tb) -> exp ( -e^tl ~ e^-tl - eEabtJb 



(4.10) 



pointwise for all |ta|,|t6| < K. Since K was arbitrary, Xn(taitb) converge every- 
where. This proves the theorem. □ 



We now finish the proof by proving Lemma 4.2 



Proof of Lemma 4-^- It is enough to prove (f.4) since the other statements follow 
immediately with Markov's inequality and Slutsky's theorem. 

We have to distinguish the case > 1 and < 9 < 1. We begin with 9 > 1 and use 
( |3.3D and the conditions (||) and (^) to get 



E 



< 



f9{9 + l) 



^(a„ + ^6™)(l;n-C(,")) 



^(|a™| + |6,„|) =0(1). 



(4.11) 



This proves the lemma for 6' > 1. 



The case < 6' < 1 is a little bit more difficult since (3.3) is now weaker. We solve 
this problem by splitting the sum. We do this as follows 



E 



Ar^ 



iBr, 



< E 



-E 



n<n/2 



(4.12) 



n/2<rn<7i 



-E 



(flm + ibm)C„ 

n/2<m<~n 

We now show that each summand in ( J:.12| ) is 0(1). 



Equation gives us E CT^ - < < ^^i^ for 1 < to < n/2. We 

thus can use the same computation as in (4.11) to see that the first summand in 
( [4.12D is equal to 0(1). 

We next look at the second summand in ( [1.121 ). We have 



E 



J2 

n/2<m<n 



< 



n/2<m<n 



\brn\ 
TO 



(4.13) 
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Wc now use Lemma 3.5 and eondition (|^) and (H) to sec that 

y \a^\ + \b„,\ ^ ^ ^ nun{5^M)_ (4.14) 



n/2<m<n 



This shows that the seeond summand in ( 4.12 ) is also 0(1). 
We finally look at the third summand. 

It is obvious that a permutation a Cz Sn can have at most one cycle with length 
greater than n/2. This fact and condition (|l|) together gives us 



E 



n /2<7n<.n 



= \a,n+ib,n\^[Cm = l] 

7i/2<7n<n 

n/2<m<n 

< J2 ^ [Cm ^ I] + 27: ^ = 1] 

7i/2<m<n n/2<m<n 



= Y \am\P[Cm^l]+2TTP 

n /2<i7n<n 

<2n+ Y = 



Y Cm>0 

n/2<m<7i 

(4.15) 



n/2<m<n 



We have used on the third line that {Cmi = 1} H {Cm2 = 1} = for mi, m2 > n/2 
and mi ^ m2. 



If the sequence {arn)m=i is bounded by a constant A', we can argue as in ( |4.15[ ) to 
see 



E 



n/2<m<n 



< 27r + is: = 0(1). 



(4.16) 



If the sequence (am,)m=i is unbounded, we have to be more careful. We first look 
at P [Cm = 1] for m > n/2. It follows immediately from (2^) that 



(7(») = i] = E Tc*"' 



\ n—m J 



^—y\ 
n } 



with 7=1 



(4.17) 



We now need an upper bound for • ^ simple computation shows 

n(n — 1) • • • (n — m + 1) 



/ n— m— 7\ 
V n—ni J 



("D (n-7)(ri-7-l)---(n-m-7 + l)' 



Wc thus have 



log 



/n—7n~j\ \ n n n 

\ n ) J fe=n— m+1 /c=n — m+1 fe=ri— m+1 

= t {h0(i)],0m., t i (4.18) 



/c— n— m+1 



fc— n— m+1 
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Wc now use J2m=i m ^ log(^i) + +0 (i) (see Theorem 3.2]) and distinguish 
the cases m — n and m < n. If m = n then we get immediately 

log(^j^]<llogin)+K2. (4.19) 



If m < n then 



/n— m— 7\ 
n— m 
7n^^ 



log I ^ "1'^, ^ I = 0(l)+7log(n)-7log(n-m) + 0( i+ ^ 



n n ^ rn 



= -7log(l--j +0(1) (4.20) 
since n — m > 0. We put everything together and get 



m — j\ 
n—m ) 



< 



K^{l-^)-\ form<n, 
K/^nP' , for m = n. 



(4.21) 



We use condition (|l|) and the Holder inequality (see Lemma 3.C) for some p, g > 1, 
specified in a moment. 



'[C„. = l] = |a„|F[C„ = l] 



E 

n/2<m<n 




'[a„ = i] 



<2A:, I- >• |a„JM I- ^ il-;rj I +^^^) ^^-^^^ 

n/2<m<n 

The second factor is a Riemann sum for - t)-^" dt. If we choose a g > 1 

with 75 < 1 then the integral exists and one can use Theorem 3.13 to see that the 
second factor converge to this integral. We now check if we can choose q in a such 
a way that the product is 0(l). We have 



79 < 1 



{i-e)< 



1 



{i-e)<i 



1 



> 



1 



p > 



1 



/2<rn<n I""*!' 



IS 



Condition (|) now ensures the existence of a p > | such that ^ J2, 
0(1). We get with this p (and q) that the product is bounded. This shows that the 
third su mma nd in ( 4.12 ) is 0(1). This prove the lemma and completes the proof of 
Theorem [O. □ 



Many assumptions we need in the proof of Theorem 4.1 are to handle the case d < I. 
If one is only interested in the case uniform measure (0 = 1) or in > 1, one can 
weaken the assumptions. We state this as a corollary. 

Corollary 4.2.1. Let 6 > I be fixed. Let (cm)m=i sequence of complex numbers 
with am = Re(c„i), b„i = Im(c,„) and 



(!') \bm\<2Tr, 

(4) ^E™=ia^->K, ^E™=i^™^^&. ^E™=iam&m 

(5) ^El^iW^\'=o{loe/'in)). 



Eab for oo, 
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Let An and be as in Theorem ^.1. We then have 



An + iBn - Ee [An + iBn] d , r 
TV 



with Af a complex normal distributed random variable with covariance matrix 

Va Eab 
Eab Vb 

The real and the imaginary part of M are independent if and only if Eab — 0. 



(4.23) 



(4.24) 



5. Proof of the main Theorem 2.7 



Wc now are ready to prove Theorem 2.7. We recommend to read first Section 3.3 
before reading tliis proof. 



Proof of Theorem 2 



Wc have by definition 

n 



(5.1) 



{mt} ,t = (im)^^i and define 



We thus can apply Theorem 4T with c,„ :— log(/(a;™)). We now show that the 
conditions (|l|) - (||) are fulfilled in all cases mentioned in Theorem 

We use the notation x = e'^^'^^,t,n 

a(s):=log|/(e2--)|, 6(s) 

/i(s) |log|/(e2--)||, k{s) 

Case 1.1 : x not a root of unity and / zero free. 

Condition (|l|) is trivially fulfilled since a™ is bounded in this case and bm < 27r by 
definition of w"{f). We next look at condition 



:arg(/(e2--)), 
Harg(/(e2"^))| 



(5.2) 
(5.3) 



-E 



E 



Re(cn 



We have 
1 " I 1 " 

-E log|/(^")| - -E'^(*" 



m— 1 



iri—X 



(5.4) 



The function /i(s) is in this case a real analytic function on [0, 1] and we therefore 
can apply Theorem 3.8 to see that the last expression converge to /i(s) ds. 

We need in d|) also the rate of convergence. We thus use Theorem [3.12 instead 
of Theorem |3.8| and have therefore to estimate -D^(t). Since x is assumed to be of 
finite type, see Definition 3.15, we can use Theorem 3.17 and get ^^(t) — 0{n~°') 
for some a > 0. This gives the desired error rate. 

It follows with the same argument that 



1 " 

— \am\^ / h^is) ds for each 1 < p < oo. 
n ^-^ In 

m=l 



(5.5) 



This shows that conditions (||), (^), (^ and the first part of condition ^ are 
fulfilled. 

We next look at bm ~ Im (log(/(a;™))) ~ arg(/(x'")) = 6(trn)- Since / is real 
analytic, one can use function theory to show that there exists a finite set D C [0, 1] 
such that b{s) is real analytic in [0, 1] \ -D and the limits linis-t-so arg(/(e^'^*'*)) and 
lims^so arg(/(e^'^*^)) exists for all sq G [0, 1]. We omit here the details since this is a 
standard argu ment. This shows that b{s) is of bounded variation and that we can 
apply Theorem 3.12 to b{s). We get 

„1 n „i 

\bm]= I \b{s)\ ds + Oin"^") and - V ^ / (^b{s)Y ds. (5.6) 



1 
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Similarly we get 



1 " 



a{s)b{s)ds= / Re(log(/(e2"^)))lm(log(/(e2--))) ds 



(5.7) 



This gives the desired expression for the covariance mentioned in Theorem | ^ 

This shows that condition an d the rest of condition (||) arc fulfilled. This com- 
pletes the proof of Theorem ^.7| in this case. 



Case 1.2 x not a root of unity and all zeros of / are roots of unity. 

The function h{s) is in this case not anymore of bounded variation and we thus 
have to apply Theorem ^.13 



The discrepancy can be estimated as in case 1.1 as Dn(t) = 0{n 

We define si, . . . , to be th e zeros of f{e'^"'^) and choose now a, S — 6{n) such 
that the error terms in ( |3.16| ) vanishes for n oo. By assumption, all zeros of / 
are roots of unity and thus there exist a q £ N such that Sk = ^ for some pk G N. 
Since x is assumed to be of finite type, we can apply Lemma 3.16 to see that 

K 



\tn 



Sk I > — for 1 < m < n,l < k < d 

'=1-^7 - - : - - 



(5.8) 



and some jjK > 0. We choose now 5 ~ Condition (|l|) follows immediately 
with this choice of 5. 



Let So — 0, Sd+i = 1 and let V{h) be as in Theorem 3.13 . We then get 













d 




[ h{t) ds 


< 


[ h{s) ds 


+ 


/ h{s) ds 




/ h{s) ds 


'o 




Jo 




Jl-S 


k=l 


J Sk-5 



+ D*^{t)V{h) + SY,\Hsk + S)\ + \h{sk+i-5) 

k=Q 

<0 (n-'^ log(n'^)) + O (n-" log(n'')) 
=0(71^''") for some Sa > 0. 

Thus condition (||) is fulfilled. A simple calculation shows that 



(5.9) 



hP{s) ^ Kp logP(s), ^h^is) ^ pKp^^^L^ 
ds s 



for 1 < p < oo. 



It is now easy to see that one can use the same argumentation as in (5.£) also for 
h?,h^ and . The conditions (|^), (||) and the first part of condition (^). The 
argumentation for bm is as above. 

Case 2 : x a root of unity of order p and /(x™) 7^ for all 1 < to < p. 
We have 



- I"™! E '^(*") = Y.Y. — pfc'^(w) 



m— 1 



Y.^^^'^) I log(n) jp + k 



(5.10) 



k=l 
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It is easy to see that 



lim 



1 1 



lim 



E 



1 1 



n->oo log(7i) ^ + A; n->oo \og{n/p) + log(p) ^ jp p 

This shows the desired convergence. To get the rate of convergence, one has to use 



V - = log(n) + Iu + 0{- 



i=i 



(5.11) 



This is a classical result and can be found for instance in the book |T|. This shows 
that Condition (|^) is fulfilled. The other calculations are similar. We therefore omit 
them. 



We have until now proven that 



-EK(/)] ^^^^ 
Vlog(n) 



(5.12) 



in all cases mentioned in Theorem p.7| inclusive the calculation of the correlation. 
To complete the proof, we have to show that 



^log(n) 



We define as in the proof of Theorem |4.1| 

n 

[J"(/)(x) := 5]log(/(x'"))y„ 



w 



(5.13) 



(5.14) 



We know from Lemma that w^{f) and uj"{f ) have the same asymptotic behav- 
ior. It is therefore enough to prove 



E[z^"(/)] 



Vlog(n)m(/) -> 0. 



(5.15) 



We first look at the case x not a root of unity. We get with Lemma |3.5| and (|5.9| ) 

E[w"{f){x)] = e ^ '^^'^ " = eiog{n) / log(/(e2--)) ds + 0{l). (5.16) 



(5.17) 



We have by definition that m(/) = log(/(e^'^'*)) ds and thus 

1 



1 I^M^ 

Vlog(n) 



In a; is a root of unity, one has to replace ( |5.9D by (|5.1l| ). We omit the details since 
this calculations are similarly. □ 



6. Estimation of the Wasserstein distance 



We estimate in this section the convergence rate of the random variable w"(/)(x) 
with respect to the Wasserstein distance, see Theorem 3.2 and Theorem 3.5. Unfor- 
tunately we have to distinguish between complex and real valued random variables. 
We thus look in Section first at the real and the imaginary part of w^{f) sepa- 
rately and then look in Section 6.2 at the complex case. 
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6.1. The real and the imaginary part separately. 

Definition 6.1. Let Xi,X2 be real valued random variables not necessarily defined 
on the same space. The Wasserstein distance dw between Xi and X2 is then defined 
as 



dw{XuX2) := sup |E [g{Xi)] - E [g{X2) 



(6.1) 



with 



(6.2) 



g := L e C^R,R) : sup \g'{t)\ < l] . 
It is easy to see that dw is a metric. We now show 

Theorem 6.2. Let 6 > 0, and x be either not a root of unity and of finite type or 
a root of unity. Let M ~ Ma + iMb be as in Theorem 2.7, then 

dw (^a, Re (^^^^^= - 0^]^)m{f){x)yj = 0(log-i/2(n)), (6.3) 
dw Ufb, Im f^^J!^= - ev/biRm(/)(x) jj = 0(log-i/2(n)). (6.4) 



We prove this theorem by reducing it to 

Theorem 6.3 (j^, Theorems 3.1 and 3.2]). Let ^1, ■ ■ ■ be independent random 
variables with 



E [U] = 0, E 



VandY,^ < 00. 



Then 



(6.5) 



dw E -^(0'^) ^ 3^3/2 ^ E [lei] 



\m—l 



711 — 1 



(6.6) 



Proof of Theorem 6^ . We use the notation 

Cm = log(/(a;™)), = Re(cm), = Re(c™), 



a 



(t)==log|/(e2"*)|, 6(0 = arg(/(e2"*)) 



(6.7) 
(6.8) 



We use as in the proof of Theorem |2/7| the Feller coupling, sec Section The 
random variables w"'{f){x) and w"'{f){x) are thus defined on the same space. We 
now get with Lemma 



dw Re 



f ^"(/)(^) ^ 



,Re 



sup 

see 



sup 

see 



E 



9 Re 



g Re 



^ ^"(/)(^) y 

'u;"(/)(x)V 



< 



«;"(/) (x) iJ"(/)(x) 



- E .9 I Re 
9 I Re 



^"(/)(a;) 
w"{f){x)\ 



Similarly for the imaginary part. 



Vlog(") 
0(log-i/2(n)). 



(6.9) 
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This shows that wc can replace w^{f) by in Theorem 3.S. It fohows with 

the triangle inequality and (5.17) that one can replace also 



We now have 

Re U"(/) - 



^yJ\og{n)m{f){x) by 



m / 



E 



(6.10) 



1 



Since all are independent, we can apply Theorem |C.3| for 

6, ' 



^Zaml^"™--). (6.11) 



.12) 



We only have to check that the assumptions of Theorem 3.3 are fulfilled. We use 
the computations in Section |^ and get 



E 



He, 



log(n) 



m— 1 

1 



2 

^^ = K + 0(n-*"), 
^ — ^ m 



m— 1 



log'/'(") i^i ™ 



ft I 



= o 



(6.13) 



(6.14) 



An applicati on of the triangle inequality shows that one can neglect the error term 
0{n~^'^) in ( |6.13 ). This proves the theorem. □ 

6.2. The two dimensional case. The definition of the Wasserstein distance can 
be extended to without any problems. Unfortunately it is often very difficult 
to handle the case c? > 1. In many situations it is much easier to take stronger 
assumptions on the test functions g. We thus set 

Definition 6.4. Let Xi,X2 be random variables with values in M.^, not necessarily 
defined on the same space. The weak Wasserstein distance d^w between Xi and 
X2 is then defined as 



with 



A^w{Xi,X2) := sup |E [g{X^)\ - E [g{X2)\\ 
g&Q 



;= {(? e C°^{M.'^,M.)-Mi{g) < 1 and Ahig) < l} 

Qk 



Mk{g) 



sup 

l<2i ,ik'^c 



• • • du^, 



(6.15) 
(6.16) 
(6.17) 



We now show 

Theorem 6.5. Let 9 > , and x be either not a root of unity and of finite type or 
a root of unity. Let M be as in Theorem 2. 1 . We then have 



?ybiRm(/)(x) =0(log-i/2(n)). 



(6.18) 



This theorem is an alternative proof of Theorem 2.7 since the weak Wasserstein 
distance is a metric on the space of random variables. 

We prove Theorem |6 . 5| with Stein's method. We can not give here a full introduction, 
but try to illustrate at least the idea of Stein's method. Assume that a random 
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variable Z and sl g Cz G arc given. In many situations one can find a "good" 
function ]j only depending on g and Z such that 

E[g{X)]-E[g{Z)]=E[giX)]. (6.19) 

for all random variables X. This simplifies the study of the Wasserstein distance 
since one now has to consider only one random variable. This reformulation is of 
course only useful if we can find a 5 with good properties. This is surprisingly 
often the case. For d = 1, a standard normal distributed random variable and 
g e CHM,K), one has 

E [f'giX) ~ Xf.iX)] = E [g{X)] - E [giN)] (6.20) 

with 

fg{x) e-'/2 r (git) - E [N])e-''^^ dt. (6.21) 



To proof Theorem 3.4, we need of course a 2-dimensional version. We use here 

Lemma 6.6 (||l2|, Lemma 1]). Let g : — )• R 6e a smooth function and M he a 
bivariate normal distributed random variable with covariance matrix E. We set 

Ugix) := J - (e [g{Vtx + VT^^e)] - E [.9(A^s)]) dt. (6.22) 

We then have 

E [giX)] - E [g{Af)] = E [{X, V(C/.g)(X)) - (Hess(Ug)(X), E)^.., ] (6.23) 
with (•, •) the standard inner product on and {Mi, M2) ^ g = Tr(MiA/J). 



Before we proof Theorem 6.5, we have to take a look at the derivations of Ug. 
Lemm a 6. 7 ( |l^ , lemma 2]). Let g : ^ M. be a smooth function and Ug as in 



Lemma 



6. t. If T, is positive definite, then 

M3{Ug)<K-M2ig) (6.24) 

with K only depending on S. 

There exists also upper bounds for Mk(Ug) if E is non negative definite, but we do 
not need them here. The reason is the following lemma. 



Lemma 6.8. The covariance matrix E in Theorem 2.1 is singular if and only if 
there exist 7 = (7a, jb) G R \ {0} such that 7a log |/(a:™)| = 76 arg (/(.t™)) for all 
m € N except finitely many. 

Proof. We prove this lemma only for x not a root of unity and of finite type. The 
case X a root of unity is similarly. We define a(s) and 6(s) as ( |6.8[ ). We know from 
Section || that 



K = / a\s) ds, Eab = I a{s)b{s) ds, H = / 
Jo Jo Jo 



{s)b{s) ds, Vb = / 6^(s) ds. (6.25) 
Jo 

Since E is a 2 x 2 matrix, one can directly compute the eigenvalues. One gets after 
a small calculation that E is non-negative definite if and only if 

a^{s) ds^ ( f b^{s) dt^>(f a{s)b{s) dt^ . (6.26) 



and E is singular if and only if we have equality in (B.26). But equation (6.26) is the 
Sehwarz inequality for L^. This shows that E is always non negative definite and 
that E is singular if and only if the functions a(t) and b{t) are linearly dependent. 
This proves the lemma since a{t) and b{t) have only finitely many discontinuity 
points and x"^^ ^ x"^^ for m\ ^ m2. □ 
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Proof of Theorem G.t. Let .g G be given with Mi{g) < 1, M2{g) < 1. We have to 
distinguish the cases S singular and S regular. 



We start with the singular case. We know from Lemma 6.8 that S is singular if and 
only if a,„ b„i for only finitely many m. We thus have 



VvlogWrt^l A 



= E 



,viogR™=i 



0(log-i/2(n)) 



with g{t) := . This shows that we can argue as in the one-dimensional case. 



We now come to S regular. 

One can use the same argumentation as in the proof of Lemma 3.2 to see that it is 
enough to show 



du,w ( A/g, , Cm [ Ym 



1 

m 



(6.27) 



with A/g a bivariate normal distribution with covariance matrix 



Va Eab 
Eab Vb 



m— 1 m Z^m— 1 jni 

En ambm sr^n b7„ 
m— 1 m Z^m— 1 m 



(6.28) 



and 



- 7 



and 



Vlog(n) 



.29) 



Wc now use ( 6.23| ) to give the desired upper bound. Let Ug he as in ( 6.22 ). We use 
the notation 



V(C/5) = ( ) , Hess(C/.9) = ( f " \ and X := V c„, f y,„ - - 



"Jab 9bb 



rn— 1 



(6.30) 



We also introduce 



m=l ^ ^ m=l ^ 



l<fc<n 
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Wc now identify C with via a + ib 



We first look at the summand 



E[{X,\/{Ug){X))]=E[XagaiX)]+E[XbgbiX)]. We have 



E[XagaiX)] = ^5„E 

m— 1 

71 



m=l 



1 

m 
1 

m 
1 

771 



fn=l 
n 

^ amE (s/ga{X^ + tc^), (j"^ ^ K 



1 

77t 



9a{X) 

(3a(X)-5a(^rn)) 

'^9a{Xrn + tCm) 



dt 



l{0<t<Ym} dt 



E 



{Xm+tCm))KrrAt) dt 

(6.31 



with 



= E 



Y 



l{o<t<y„j 



.32) 



We have used for the second equahty that E \Y„i — 1 /777] = and that X^ is inde- 
pendent of Ym- We have of course also to justify the existence of the integrals, but 
this follows immediately from Mi{g) < l,M2{g) < 1 and Lemma 3.7. 
We next look at 

E [(Hess(C/.g)(X), S)^ = KE [gaa{X)] + 2EabE [gabiX)] + HE [gbbiX)] . (6.33) 
A direct computation shows that 

/ K^{t) dt = — and / tKm{t) dt = — + O ( ^ ) (6.34) 
Jo ™ Jo V"^ / 

with O(-) independent of n and S. We thus have 



E 



VagaaiX) 



= E 



■JaaiX) V al^ — 



rn— 1 



5™ / gaa{X)K„,{t) dt 



(6.35) 



We combine ( |6.31 ) and (6.35) and get 

E [(X, V(;7.9)(X))] - E [(Hess([/g)(X), I])^.^.] 



am{9aa{Xm + Cmt) - gaa{Xm + Cm 

+ ^ 2E / {Xm + Cmt) ~ gab{Xm + C„iYm) ) A; 

m=l '-"'0 
n p „ 



.(t) dt 



6„ gbb{Xm + Cmt) - gbb{Xm + CmYm) Km{t) dt 



.36) 
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Wc now use ( ^.34 ) and Lemma 6.7 to get 
E 



(^gaa{Xm + Cmt) - gaa{Xm + CmYm)^Km{t) dt 

M3{Ug)\d„,\{t + Y„,)K„,{t) dt 



< 



< M2{g)\c^\ / {t+-)K„M dt 
m 



< 



K 



-.AUg)'- 



I ~ I ^771 1 



Thus 



Vlog(n) 

; [{X, V(C/.g)(X))] - E [(Hess(C/g)(X), 1])^,^.; 



(6.37) 
(6.38) 



< 



K 



It follows with a computation similar to computation in ( |5.9| ) that the last expression 
is 0(log~^/^(n)) . This proves the theorem. 



□ 
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